Elliptic solutions to difference non-linear equations and 
nested Bethe ansatz equations 

I.Krichever * 
March 2, 1998 



Abstract 



We outline an approach to a theory of various generalizations of the elliptic Calogero- 
Moser (CM) and Ruijsenaars-Shneider (RS) systems based on a special inverse problem 
for linear operators with elliptic coefficients. Hamiltonian theory of such systems is 
developed with the help of the universal symplectic structure proposed by D.H. Phong 
and the author. Canonically conjugated action-angle variables for spin generalizations 
of the elliptic CM and RS systems are found. 
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1 Introduction 



The elliptic nested Bethe ansatz equations are a system of algebraic equations 

a(x? - x] +1 )a(x? - 77 - aff )g(s? - x]- 1 + 77) _ i 
W - x7 _1 )t7(x? + 77 - - a;? +1 - rf) 1 ' J 

for iV unknown functions Xj = xf, i = 1, . . . , N, of a discrete variable n (Q). (Here and 
below cr(x) = cr(x\ui,U2), ((x) = a/), and p(x) = p(x\u, a/) are the Weierstrass a-, 

and p-functions corresponding to the elliptic curve with periods 2a; , 2a;'.) This system is an 
example of the whole family of integrable systems which have attracted renewed interest for 
years. The most recent burst of interest is due to the unexpected connections of these systems 
to Seiberg-Witten solution of N = 2 supersymmetric gauge theories @. It turns out that 
the low energy effective theory for SU(N) model with matter in the adjoint representation 
(identified first in || with SU(N) Hitchin system) is isomorphic to the elliptic CM system. 
Using this connection quantum order parameters were found in [||]. 

The elliptic Calogero-Moser (CM) system [|j], || is a system of N identical particles on a 
line interacting with each other via the potential V(x) = p(x). Its equations of motion have 
the form 

x i = AY^p\x i -x j ). (1.2) 

The CM system is a completely integrable Hamiltonian system, i.e. it has N independent 
integrals in involution (0, ||). The second integral H 2 is the Hamiltonian of ( |1.2|) . 

In a remarkable connection of the CM system with a theory of elliptic solutions to the 
KdV equation was revealed. It was shown that the elliptic solutions of the KdV equations 
have the form u(x,t) = 2X4=1 p( x ~ x i(t)) an d the poles Xi(t) of the solutions satisfy the 
constraint J2j^i p'( x i ~ x j) = 0, which is the locus of the stationary points of the CM system. 
Moreover, it turns out that the dependence of the poles with respect to t coincides with 
the Hamiltonian flow corresponding to the third integral H 3 of the system. In |[LU|| , [IT| it 



was found that this connection becomes an isomorphism in the case of the elliptic solutions 
to the Kadomtsev-Petviashvilii equation. Since then, the theory of the CM system and its 
various generalizations is inseparable from the theory of the elliptic solutions to the soliton 
equations. 



In [0 system (|1.1|) was revealed as the pole system corresponding to the elliptic solutions 
of completely discretized version of the KP equation on lattice. It was noticed that equations 
J]) have the form of the Bethe ansatz equations for the spin-^ Heisenberg chain with 



impurities. Its connection with the nested ansatz equations for the A^-lattice models was 
established in M. 



As shown in [0 an intermediate discretization of the KP equation which is the 2D Toda 



lattice equations leads to the Ruijesenaars-Schneider system |I3| : 



' x i = ^2xiX s (V(xi - x s ) - V(x s - Xi)), V(x) = ((x) - C(x + 77), (1.3) 



which is a relativistic version of (|1.2|) 



2 



The main goal of this paper is to present a general approach to the theory of the CM type 
many body systems which is based on a special inverse problem for linear operators with the 
elliptic coefficients. This approach originated in fTS] and developed in M, flT5[, [ITB] clarifies 



the connection of these systems with the soliton equations for which the corresponding 
linear operator is the Lax operator. We formulate the inverse problem in the next section 
and show that its solution is equivalent to a finite-dimensional integrable system. We discuss 
an algebraic-geometric interpretation of the corresponding systems. 

The advantage of our approach is that it generates the finite-dimensional system simulta- 
neously with its Lax representation. Until recently, among its disadvantages was the missing 



connection to the Hamiltonian theory. For example, in spin generalization of the RS 
system was proposed and explicitly solved in terms of the Riemann theta-functions of aux- 
iliary spectral curves. At the same time, all direct attempts to show that this system is 
Hamiltonian have failed, so far. 

One of the existing general approaches to the Hamiltonian theory of the CM type systems 
is based on their geometric interpretation as reductions of geodesic flows on symmetric spaces 
f|. Equivalently, these models can be obtained from free dynamics on a larger phase space 



possessing a rich symmetry by means of the Hamiltonian reduction JT7J. A generalization 
to infinite-dimensional phase spaces (cotangent bundles to current algebras and groups) was 
suggested in f[8 |, [JT9| . The infinite-dimensional gauge symmetry allows one to make the 



reduction to a finite number of degrees of freedom. 

A further generalization of this approach consists in considering dynamical systems on 
cotangent bundles to moduli spaces of stable holomorphic vector bundles on Riemann sur- 
faces. Such systems were introduced by Hitchin in the paper fl20 |, where their integrability 



was proved. An attempt to identify the known many body integrable systems in terms of 



the abstract formalism developed by Hitchin was made in [pi]] . To do this, it is necessary to 
consider vector bundles on algebraic curves with singular points. It turns out that the class 
of integrable systems corresponding to the Riemann sphere with marked points includes spin 



generalizations of the CM model as well as integrable Gaudin magnets (see also [23|). 

Unfortunately, a geometric interpretation of spin generalization of the elliptic RS system 
has not been yet found. Recently, such realization and consequently the Hamiltonian theory 
of the rational degeneration of that system were found in |24| . 



In section 3 we develop Hamiltonian theory of the CM type systems in the framework of 
the new approach to the Hamiltonian theory of soliton equations proposed in and ||26|| . 
It can be applied evenly to any equation having the Lax representation. The symplectic 
structure is constructed in terms of the Lax operator, only. We discuss three basic examples: 
spin generalizations of the CM and the RS systems, and the nested Bethe ansatz equations. 
We would like to emphasize that the universal form of the symplectic structure provides a 
universal and direct way to the action-angle type variables. 

We would like to refer to || for the analysis of connections of the elliptic CM system to 
Seiberg-Witten theory of N = 2 supersymmetric gauge theories. 
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2 Generating problem 



Let £ be a linear differential or difference operator in two variables x, t with coefficients which 
are scalar or matrix elliptic functions of the variable x (i.e. meromorphic double-periodic 
functions with the periods 2u a , a = 1, 2). We do not assume any special dependence of the 
coefficients with respect to the second variable. Then it is natural to introduce a notion of 
double-Bloch solutions of the equation 

C^ = 0. (2.1) 

We call a meromorphic vector-function f(x) that satisfies the following monodromy proper- 
ties: 

f{x + 2u a ) = B a f{x), a = 1,2, (2.2) 

a double-Bloch function. The complex numbers B a are called Bloch multipliers. (In other 
words, / is a meromorphic section of a vector bundle over the elliptic curve.) 

In the most general form a problem that we are going to address is to classify and to 



construct all the operators C such that equation ( |2.1| ) has sufficiently enough double-Bloch 
solutions. 

It turns out that existence of the double-Bloch solutions is so restrictive that only in 
exceptional cases such solutions do exist. A simple and general explanation of that is due to 
the Riemann-Roch theorem. Let D be a set of points X{, i — 1, . . . , m, on the elliptic curve 
T with multiplicities di and let V = V(D; Bi, B 2 ) be a linear space of the double-Bloch 
functions with the Bloch multipliers B a that have poles at Xi of order less or equal to di and 
holomorphic outside D. Then the dimension of D is equal to: 

dim D = deg D = dj. 

i 

Now let Xi depend on the variable t. Then for / G D(t) the function Lf is a double-Bloch 
function with the same Bloch multipliers but in general with higher orders of poles because 
taking derivatives and multiplication by the elliptic coefficients increase orders. Therefore, 
the operator £ defines a linear operator 

C\ D : V(D(t);B 1 ,B 2 ) i— > V(D'(t); B u B 2 ), N' = deg/T > N = deg A 

and (|2.1| ) is always equivalent to an over- determined linear system of N' equations for iV 
unknown variables which are the coefficients c, = q(t) of an expansion of ^ G V(t) with 
respect to a basis of functions fi(t) G V(t). With some exaggeration one may say that in the 
soliton theory the representation of a system in the form of the compatibility condition of an 
over-determined system of the linear problems is considered as equivalent to integrability. 

In all of the examples which we are going to discuss N' = 2N and the over-determined 
system of equations has the form 

LC = kC, d t C = MC, (2.3) 
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where L and M are N x N matrix functions depending on a point z of the elliptic curve as 
on a parameter. A compatibility condition of ( [2.3|) has the standard Lax form d t L = [M, L], 
and is equivalent to a finite-dimensional integrable system. 

The basis in the space of the double-Bloch functions can be written in terms of the 
fundamental function <$>(x,z) defined by the formula 

= ^l=^e«*. (2.4) 
a[z)a{x) 

Note, that z) is a solution of the Lame equation: 

(J^ - 2p(x)j ${x, z) = p(z)$(x, z). (2.5) 

From the monodromy properties it follows that $ considered as a function of z is double- 
periodic: 

$(x, z + 2cj q ) = $(x, z), 

though it is not elliptic in the classical sense due to an essential singularity at z = for 
x ^ 0. 

As a function of x the function <&(x, z) is double-Bloch function, i.e. 

$(x + 2u a , z) = T a (z)$(x, z), T a (z) = exp (2u a ((z) - 2((u a )z) . 

In the fundamental domain of the lattice defined by 2u a the function $(x, z) has a unique 
pole at the point x = 0: 

= a;" 1 + 0(x). (2.6) 

The gauge transformation 

f{x) ' ^ /(x) = /(x)e M , 

where a is an arbitrary constant does not change poles of any function and transform a 
double Bloch-function into a double-Bloch function. If B a are Bloch multipliers for / than 
the Bloch multipliers for / are equal to 

B 1 = B x e 2au) \ B 2 = B 2 e 2aw2 . (2.7) 

The two pairs of Bloch multipliers that are connected with each other through the relation 
( j2.7p for some a are called equivalent. Note that for all equivalent pairs of Bloch multipliers 
the product B^B^ 1 is a constant depending on the equivalence class, only. 

From ( p.6|) it follows that a double-Bloch function f(x) with simple poles Xi in the 
fundamental domain and with Bloch multipliers B a (such that at least one of them is not 
equal to 1) may be represented in the form: 

N 

f(x)=J2c i Hx-x l ,z)e kx , (2.8) 
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where Cj is a residue of / at Xj and z, k are parameters related by 

B a = T a (z)e 2u) « k . (2.9) 



(Any pair of Bloch multipliers may be represented in the form (|2.9| ) with an appropriate 
choice of the parameters z and k.) 

To prove (|2.8| ) it is enough to note that as a function of x the difference of the left and 
right hand sides is holomorphic in the fundamental domain. It is a double-Bloch function 
with the same Bloch multipliers as the function /. But a non-trivial double-Bloch function 
with at least one of the Bloch multipliers that is not equal to 1, has at least one pole in the 
fundamental domain. 

Now we are in a position to present a few examples of the generating problem. 



Example 1. The elliptic CM system ([[L4 
Let us consider the equation 

C^j = (d x -dl + u{x,t))y = 0, (2.10) 

where u(x,t) is an elliptic function. Then as shown in [|lj|] equation (|2.10| ) has A" linear 
independent double-Bloch solutions with equivalent Bloch multipliers and N simple poles at 
points Xi(t) if and only if u(x,t) has the form 

N 

u{x,t) = 2j2p(x-Xi{t)) (2.11) 

i=l 

and Xi(t) satisfy the equations of motion of the elliptic CM system ( p..2| ). 

The assumption that there exist A" linear independent double-Bloch solutions with equiv- 
alent Bloch multipliers implies that they can be written in the form 

JV 

* = £ k, z)Q(x - Xi {t), z)e kx+kH , (2.12) 

i=l 

with the same z but different values of the parameter k. 



Let us substitute ( |2.12j ) into ( |2.10| ). Then ( |2.10| ) is satisfied if and if we get a function 



holomorphic in the fundamental domain. First of all, we conclude that u has poles at x,, 
only. The vanishing of the triple poles (x — Xj)~ 3 implies that u(x,t) has the form ( 2.11|) . 



The vanishing of the double poles (x — x,i) 2 gives the equalities that can be written as a 
matrix equation for the vector C = (cj): 

(L(t,z) -kI)C = 0, (2.13) 

where / is the unit matrix and the Lax matrix L(t,z) is defined as follows [] : 

L i:j (t, z) = --SijXi - (1 - 6ij)$(xi - Xj, z). (2.14) 



bn order to simplify the consequent formulae in Section 3 and present all the examples in the same 
framework we choose here and in the next example a normalization of L which differs from that used in Jl4| , 
H by the factor -1/2. 
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Finally, the vanishing of the simple poles gives the equations 



(d t -M(t,z))C = 0, (2.15) 

where 

Mij = (p(z) -2J2 P(*i ~ Xj) j 6ij - 2(1 - <%) $'0; - Xj, z). (2.16) 

The existence of A" linear independent solutions for (|2.10|) with equivalent Bloch multipliers 
implies that ( |2.13|) and ( p,15|) have N independent solutions corresponding to different values 
of k. Hence, as a compatibility condition we get the Lax equation L = [M, L] which is 
equivalent to (|L~2|) . Note that the last system does not depend on z. Therefore, if Q2.13|) 
and ( |2.15| ) are compatible for some z then they are compatible for all z. As a result we 
conclude that if ( p.lOj ) has A" linear independent double-Bloch solutions with equivalent 
Bloch multipliers then it has infinitely many of them. All the double-Bloch solutions are 
parameterized by points of an algebraic curve T defined by the characteristic equation 

N 

R(k, z) = det{kl - L{z)) = k N + ri{z)k N -* = 0. (2.17) 

i=l 

Equation ( |2.17[ ) can be seen as a dispersion relation between two Bloch multipliers and 
defines T as A^-sheet cover of IV 

From (U) and ([HJ) it follows that 

L = GLG-\ Gij = e az)xi Sij. (2.18) 

At z = we have the form 

L = z-\F-I) + 0(1), Fy = l. (2.19) 

The matrix F has zero eigenvalue with multiplicity A^ — 1 and a simple eigenvalue A^. 
Therefore, in a neighborhood of z = the characteristic polynomial (|2.17|) has the form 

N 

R(k,z) = H(k + nz- 1 +hi + 0(z)), ut = 1 - N, Vi = 1, i > 1. (2.20) 

i=l 

We call the sheet of the covering r at z = corresponding to the branch k = z~ 1 (N— 1)+0(1) 
by upper sheet and mark the point Pi on this sheet among the preimages of z = 0. From 
( f2.20| ) it follows that in general position when the curve T is smooth, its genus equals N. 

The coefficient Vi(z) in (|2 . 1 7|) is an elliptic function with a pole at z = of order i. The 
coefficients Jj )S of the expansion 

i-2 



r t (z)=I tfi + Y,ksd s Mz) (2.21) 



s=0 
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are integrals of motion. From ( |2.20| ) it is easy to show that there are only N independent 
among them. Recently, a remarkable explicit representation for the characteristic equation 
( FT7Q was found (§: 

R(k, z) = f(k - C(*), *), /(*, *) = C7 ^ + (*), (2.22) 

where H(k) is a polynomial. Note that ( j2.22j ) may be written as: 



/(M 



a 



1 * 1 f d\ n 



The coefficients of the polynomial are free parameters of the spectral curve of the CM 
system. 

When the Lax representation and the corresponding algebraic curve V are constructed, 
the next step is to consider analytical properties of the eigenvectors of the Lax operator L 

on r. 

Let L be a matrix of the form (|2.14| ). Then the components of the eigenvector C(P) = 
(q(P)), P = (k,z) G T, normalized by the condition 

N 

5>(P)$(-a;,,z) = i, (2.23) 

8=1 

are meromorphic functions on T outside the preimages Pi of z = 0. They have N poles 
7i,...,7tv (in a general position 7 S are distinct). At the points Pj the coordinates q(P) 
have the expansion: 

d = z{c\ + 0{z))e^ z) ^] a = {4 + 0(z))e^ z)x \ i > 1. (2.24) 

If we denote the diagonal elements of L by —pi/2, then the above constructed correspondence 

Pl , x t v— >{T,V = { ls }} (2.25) 

is an isomorphism (on the open set). 

Now let Xi(t) be a solution of (|1.2|) then the divisor V corresponding to Pi = Xi(t), Xi(t) 
depends on t, T> = T>(t). It turns out that under the Abel transform this dependence 
becomes linear on the Jacobian J(T) of the spectral curve. The final result is as follows. 

Theorem 2.1 The coordinates of the particles Xi(t) are roots of the equation 

9(Ux + Vt + Z ) = 0, (2.26) 

where = 9(£\B) is the Riemann theta-function corresponding to matrix of h-periods of 
holomorphic differentials on T; the vectors U and V are the vectors of b-periods of normalized 
meromorphic differentials on T, with poles of order 2 and 3 at the point Pi; the vector Z is 
the Abel transform of the divisor T> = T>(0). 
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This result can be reformulated in the following more geometric way. Let J(T) be the 
Jacobain (iV-dimensional complex torus) of a smooth genus N algebraic curve T. Abel 
transform defines imbeding of T into J(T). A point P G T defines a vector U in J(T) that is 
the tangent vector to the image of T at the point. Let us consider a class of curves having the 
following property: there exists a point on the curve such that the complex linear subspace 
generated by the corresponding vector U is compact, i.e. it is an elliptic curve r . This 
means that there exist two complex numbers 2u a , Im 002/^1 > 0, such that 2u a U belongs to 
the lattice of periods of the holomorphic differentials on T. From pure algebraic-geometrical 
point of view the problem of the description of such curves is transcendental. It turns out 
that this problem has an explicit solution, and algebraic equations that define such curves are 
the characteristic equation for the Lax operator corresponding to the CM system. Moreover, 
it turns out that in general position Tq intersects theta-divisor at iV points X{ and if we 
move Tq in the direction that is defined by the vector V of the second jet of T at P then the 
intersections of Tq with the theta-divisor move according to the CM dynamics. 

Example 2. Spin generalization of the elliptic CM system ([Q). 

Let C be an operator of the same form (|2.10|) as in the previous case, but now u = u^(x, t) 
is an elliptic (I x I) matrix function of the variable x. We slightly reformulate the results of 
|15| to a form that would be used later. 

Equation ( |2.10| ) has N > I linear independent double-Bloch solutions with N simple 



poles at points Xi(t) and such that (I x N) matrix formed by its residues at the poles has 
rank I if and only if: 

(i) the potential u has the form 

N 

u = "£a i (t)bt(t)p(x-x i (t)), (2.27) 

where = (ai. a ) are /-dimensional vectors and bf = (bf) are /-dimensional co-vectors; 

(ii) Xi{t) satisfy the equations 

*i = Y.( h ta 3 )^a i )p'{x i - Xj ); (2.28) 

(iii) the vectors a« and co-vectors bf satisfy the constraints 

1 

bt ai = Y,b?(t)a ta (t) = 2 (2.29) 

Q = l 

and the equations 

(H = - ^2 aj(bfai)p(xi - Xj) - A^i, b { = ^ b j (bfa j )p(x i - Xj) + Xfif (2.30) 

where Aj = \i(t) are scalar functions. 

In order to get these results we represent a double-Bloch vector function \1/ in the form: 

N 

V = J2 Si(t, k, z)*(x - Xi(t), z)e kx+k2 \ (2.31) 
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where are /-dimensional vectors and substitute it into Q2.10D . The vanishing of the triple 
pole (x — Xi)~ 3 implies that u has the form ( |2.27 ), the vectors Sj are proportional to a^, i.e. 
Sj = CjOj, where c$ are scalars, and that the constraints fl2.29| ) are fulfilled. 

The vanishing of the coefficients in front of (x — Xi)~ 2 implies that the vector C with the 
coordinates q satisfies ( ]2.13|) , where the Lax matrix has the form 

U 3 = -\8 ljXi - 1(1 - *«)(tfa,)*(*< - Xj , z). (2.32) 



The vanishing of the coefficients in front of (x — xf) 1 implies ( 2.30 ) for the vectors a« and 
the equation ( 2.15| ), where the matrix M has the form 

M tj = (Xi + tpizWy - (1 - SyW+aA&ixi - xj). (2.33) 

The Lax equation for these matrices implies ( |2.28| ) and equations ( |2.30| ) for bf (here we use 
the assumption that d{ span the whole /-dimensional space). 

System (|2.30|) after the gauge transformation 

Oj -> ciiqi, bf -> hq~ l , q { = exp (J \(t)dtj , (2.34) 
which does not effect (|2.28|) and Q2.29|) , becomes 

di = -J2 a>jQ>fai)p(xi - xj), k = Yl b j( b t a j)p( x i ~ x j)- ( 2 - 35 ) 

Equations ( p. 29|) , ( p.35 ) are invariant under the transformations 

a t ^\- l a u bf^Xibf, a^W^a^ bf -+ bfW } (2.36) 

where Aj are constants and W is a constant (/ x /) matrix. A factorization with respect to 
these transformations leaves us with a reduced phase space M. of the dimension dim Ai = 
2NI — /(/ — 1). Let us introduce a canonical system of coordinates on Ai. First of all, for 
any set of dj, bf we define the matrix 

N 

Sf = £a ia fcf ( 2 - 37 ) 
i=i 

and diagonalize it with the help of the matrix Wl a which leaves the co- vector (1, . . . , 1) 
invariant, i.e. 

SiW*p = 2KjWL £ W L = 1, 3 = 1, • • • ,1 (2-38) 

a 

Then we define 

Ai = W^a*, Bf = bfW . (2.39) 
The vectors Ai and co-vectors Bi satisfy the conditions 

J2A ta B^ = 2^X, (2.40) 

i=l 
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which destroy the second half of the gauge transformations ( |2.36| ). 
At z = we have 

L = GLG-\ L = z-\F - /) + 0(1), F v = ~(fc+ % ) = i(fljMi). (2.41) 



where G is the same as in (|2.18| ). The matrix F has rank I. Its null subspace is a subspace 
of vectors such that 



N 



Ai t aCi — 0. 



(2.42) 



i=l 



Relations (|2.40|) imply that the eigenvector of F corresponding to non-zero eigenvalue 2Kj is 
identified with Bf. 



From ( |2.41|) it follows that 



N 



R(k, z) = Y\_(k + ViZ 1 + hi + O(z)), f; = 1 — «», i < i, i/j = 1, z > i. 



(2.43) 



i=l 



As shown in JT5] expansion ( 2.43|) implies that the spectral curve defined by ( [2.17|) has (in 
general position) genus g = Nl — 1(1 + l)/2 + 1. At the same time ( 2.43[ ) implies that a 
number of independent integrals given by ( |2.17|) is equal to |dim M.. 

The angle-type variables of our reduced system are the divisor of poles of the solution of 
( P- 13] ) with the following normalization: 



I N 



Y Y A ia Ci§(-Xi, z) = 1. 



a=l i=l 

At the points P~ the components of C has the form 



d = z(4 + 0{z))e« z >\ i<l;ci= (4 + 0{z))e^ x % j > I, 



where 



jV 



N 



(2.44) 

(2.45) 
(2.46) 



Y A ita 4 = 8 J a , j<l, Y A i,*4 = o, j > i 
1=1 1=1 

The last formulae identify the normalization Q2.44j) with a canonical normalization used in 
the soliton theory. The coordinates 



N 



^ ] Aj a Cj^?{x .'i'.j. Z^ 



kx 



8=1 



of the corresponding Baker-Akhiezer function \1/ are meromorphic outside the punctures 
Pj, j < I. At Pj they have the form: 



^ Q =(5i + 0(z))e z ~ 1 ^ 



(see details in 15 
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Outside the punctures Pj the canonically normalized vector C has g + I — 1 poles 
7i, . . . , 7 g _z + i. This number is equal to | dim .M. Note that these poles are indepen- 
dent on the transformations ( |2.36D . Therefore, we have constructed an algebraic-geometric 
correspondence 

M = {x t , Pl , A h B t } ^{T,V= ( 7s )} , (2.47) 

which is an isomorphism on the open set. The reconstruction formulae for solutions of 
the elliptic CM system ( |2~28l) , ( |2~29"D , (|2~35p can be found in 0- lt turns out tnat tlie 
coordinates of Xi are defined by the same equation ( |2.26|) . The only difference is a set of the 
corresponding algebraic curves that are defined by the characteristic equation for the Lax 
matrix L of the form (|2.32|) . In pure algebraic-geometric form they can be described in a 
way similar to that in the previous example. Namely, they are curves such that there exists 
a set of / points on it with the following property: a linear space spanned by the tangent 
vectors to the curve at these points in the Jacobian, contains a vector U which spans in J(T) 
an elliptic curve r . 

Example 3. Spin generalization of the elliptic RS system (|TO) 
Let us consider now the differential-difference equation 

Cm = d t $?(x, t) - ^(x + Tj,t)- v{x, t)V(x, t) = 0, (2.48) 

where rj is a complex number and v(x,t) is an elliptic (/ x I) matrix function. It has iV > I 
linear independent double-Bloch solutions with N simple poles at points Xi(t) and such that 
(/ x N) matrix formed by its residues at the poles has rank I if and only if: 

(i) the potential u has the form 

N 

v = Y / a i (t)bt(t)V(x-x i (t)), V = ((x-x i + r } )-({x-x i \ (2.49) 

8=1 

where /-dimensional vectors and bf = (bf) are /-dimensional co- vectors; 

(ii) Xi(t) satisfy the equations 

Xl = XX&^X&ia;) (y( Xi - x,) - V( Xj - x t )) ■ (2.50) 

(iii) the vectors a« = (a^) and co- vectors 6j = (bf) satisfy the constraints 

i 

bta i = J2 b t(t)^ a (t) = x l (2.51) 

a=l 

and the system of equations 

a i = H aj(bfai)V(xi - xj) - AjOj, h = - ^ bj{bf a j) v { x j ~ x i) + '> (2.52) 

where Aj = \i(t) are scalar functions. 

The gauge transformation ( p.36| ) allows us to eliminate A, in (|2.52| ). The corresponding 
system was introduced in Jl6| and is spin generalization of the elliptic RS model, which 
coincides with ( |2.50[ ), ( |2.51| ) for / = 1. 
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Construction of integrals of system ( |2.50| )-( p.52| ) and the angle type coordinates is par- 
allel to the previous case but requires some technical modification because £ is a difference 
operator in x. First of all, we choose a different basis in a space of double-Bloch functions. 
The proper choice is defined by the formula 



§(x, z) 



a(z + x + rj) 
a(z + r])a(x) 



a[z 



rj) 



-i x/2-q 



o\z 



v) 



(2.53) 



It satisfies the difference analog of the Lame equation ( |2.5|) : 



§(x + rj, z) + c(x)$(x -T),z) — E(z)$(x, z), 



where 



c(x) 



<t(x — rj)a{x + 2rj) 



E(z) 



a{2rj) 



o{z\ 



a(x + i])a(x) a(rj) (cr(z — rj)a(z + rj)) 1 / 2 

The Riemann surface T of the function E(z) is a two-fold covering of the initial elliptic 
curve T with periods 2uj a , a = 1, 2. Its genus is equal to 2. 

As a function of x the function z) is double-Bloch function. In the fundamental 
domain of the lattice defined by 2u a , the function z) has a unique pole at the point 
x = 0: 

.-1 , A , ^ A _ H- , ^ , 1 1_ a ( Z -V) 



z) 



X 



+ A + 0(x), A 



C(z + r]) + —\n^ 
2r] o(z 



Therefore, we may represent a double-Bloch solution ^/ of fl2.48|) in the form: 

N 

* = Y,^k,z)$(x-Xi(t),z)k x/ \ 



(2.54) 



(2.55) 



i=l 



substitute this ansatz into the equation and proceed as before. We get (|2.50|) - (|2.52p . The 
corresponding Lax operators have the form 



Lij(t, z) 



~aj)§(xi 



r),z) 



Mij(t, z) = (Aj - (((rj) - A)xi)Sij + (1 - Sij)^taj)^{xi - x h z). 



(2.56) 
(2.57) 



Explicit formulae in terms of the Riemann theta functions are the same as for spin general- 
ization of the CM system. The only difference is due to the different family of the spectral 
curves. In the case I = 1 they may be defined as a class of curves having the following 
property: there exists a pair of points on the curve such that the complex linear subspace 
spanned by the corresponding vector U is an elliptic curve To- If we move To in the direction 
that is defined by the vector V + (V~) tangent to T G J{T) at the point P + (P~), then the 
intersections Xi of To with the theta-divisor move according to the RS dynamics. The spec- 
tral curves for I > 1 are characterized by the existence of two sets of points P^ 1 , % — 1 , . . . , I 
such that in the linear subspace spanned by the vectors corresponding to each pair there 
exists a vector U with the same property as above. 
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Example 4. The nested Bethe ansatz equations ([[L]]) 
Let us consider two-dimensional difference equation 

$>(x, m + 1) = + 77) + v(x, m)ty(x, m) (2.58) 

with an elliptic in x coefficient of the form 

v(x, m) = [[ — x — — . 2.59 

t=i a \ x — x" )a(x — x™ + fj) 

It has N linear independent double-Bloch solutions with equivalent Bloch multipliers if and 
only if the functions x™ of the discrete variable m satisfy the nested Bethe ansatz equations. 
The Lax representation for this system has the form 

L(m + l)M(m) = M(m)L(m + 1), (2.60) 

where 

Ly(m) = Ai(m)$(af - xf - V ,z), My (m) = ^(m)$(:zf +1 - x™, z) (2.61) 

and 



nf =1 ,^ - xf) nf =1 - - ??) 



A *(™) = TTM Z^TRn _,_ m T^Zl IT » ( 2 - 62 ) 



A class of the spectral curves is the same as for the RS system. The solution xf of fll.lj ) 
corresponding to the spectral curve and the divisor on it is defined by the equation 

6(Ux + Vm + Z) = 0. (2.64) 

Here V is the vector from the puncture P + to the third point Q on V. When this point tends 
to P + the vector V becomes the tangent vector to the curve and we come to the RS system 
as a continuous limit of ( [L.l|) . 



3 Hamiltonian theory of the CM type systems. 

As we have seen various generating problems lead to various integrable finite-dimensional 
systems which can be explicitly solved via the spectral transform of a phase space Ai to 
algebraic-geometric data. On this way we do not use Hamiltonian description of the system. 
Moreover, a'priori it's not clear, why all the systems which can be constructed with the help 
of the generating scheme are Hamiltonian. In this section we clarify this problem using the 
approach to the Hamiltonian theory of soliton equations proposed in and developed in 



First of all, let us outline a framework that was presented in the previous section. The 
direct spectral transform identifies a space of solutions with a bundle over a space of the 
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corresponding spectral curves. The fiber over a curve T is a symmetric power S 9+ T (i.e. 
an unordered set of (g + 1 — 1) points 7 S 6 T). Dimension of the space of the spectral curves 
equals to g + 1 — 1 = \ M.. The spectral curves are realized as iV-sheet covering of the elliptic 
curve T . 

Let us consider the case of spin generalization of the CM system. Entries of L(z) are 
explicitly defined as functions on j\A. Therefore, L(z) can be seen as an operator- valued 
function and its external differential 5L as an operator- valued one-form on M.. Canonically 
normalized eigenfunction C(z,k) of L(z) is the vector- valued function on Ai. Hence, its 
differential is a vector-valued one-form. Let us define a two-form on A4. by the formula 



JV 

uj = 



N 

J>es Pi < C*(z,k)(6L(z) + 5k) A5C(z) > dz, (3.1] 



where C*(z,k) is the eigen-covector (row vector) of L(z), i.e. the solution of the equation 
C*L = kC*, normalized by the condition < C*(z, k)C(z, k) >= 1. The form uj can be 
rewritten as 

uj = res 2=0 Tr (p- 1 {z)5L(z) A 5C(z) - C~\z)5C{z) A 5k) dz, (3.2) 

where C(z) is a matrix with columns C(z,kj); kj = kj(z) are different eigenvalues of L(z) 
and k(z) is the diagonal matrix kj(z)5ij. 



Remark 1. The right hand side of formula ( |3.1|) is not gauge invariant. In [£5[ a gauge 
(Ci = gC, Li = gLg^ 1 ) was chosen in such a way that uj is equal to the sum of residues of 
the differential < C{5Li A 5C X > dz. 

Note that C* are rows of the matrix C^ 1 . That implies that C* as a function on the 
spectral curve is: meromorphic outside the punctures; has poles at the branching points of 
the spectral curve, and zeros at the poles 7 S of C. These analytical properties are used in 
the proof of the following theorem. 

Theorem 3.1 The two-form uj equals 

9+1-1 



UJ 



2 5z( ls )A5k( ls ). (3.3) 



The meaning of the right hand side of this formula is as follows. The spectral curve by 
definition arises with the meromorphic function k(Q) and multi- valued holomorphic function 
z(Q). Their evaluations k(j s ), z(j s ) at the points 7 S define functions on the space j\4, and 
the wedge product of their external differentials is a two-form on M.. Formula ( |3.3| ) identifies 
k s = k(j s ), z s = z(7a) as Darboux coordinates for uj. 



Remark 2. The right hand side of ( |3.3| ) can be identified with a particular case of universal 
algebraic-geometric symplectic forms proposed in |25[| . They are defined on the generalized 
Jacobian bundles over a proper subspaces of the moduli spaces of Riemann surfaces with 
punctures. In the case of families of hyperelliptic curves that form was pioneered by Novikov 



and Veselov \2( 
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Remark 3. Equations ( |2.28| ), ( [2.291 ), ( j2.35| ) are linearized by generalized Abel transform 



A : S^^T i — > J(T) x C 1 - 1 . (3.4) 

This transform is defined with the help of a basis of the normalized holomorphic differentials 
doJii % < g, and with the help of normalized meromorphic differentials dQj, j < (I — 1) of the 
third kind with residues 1 and —1 at the points Pj and Pi, respectively. Let Q be a point of 
T. Then we define (g + 1 — 1) -dimensional vector A k {Q) with the coordinates 

rQ rQ 
Ai(Q) = I du u A g+j {Q) = I dflj. 

The isomorphism (|3.4|) is given by 

9+1-1 

<j) k = J2 A k(pts)- 

The action variables I k canonically conjugated to <fi, i.e. such that 

g+i-i 
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$4>k A 5I k , 

k=l 



are equal to 

Ii = f kdz, I g+ j = respj kdz = —fj, 

Jcft 

where a° are a-cycles of the basis of cycles on T with the canonical matrix of intersections. 

Let us outline the proof of Theorem 3.1. The differential Q =< C*5L A 5C > dz is 
a meromorphic differential on the spectral curve (the essential singularities of the factors 
cancel each other at the punctures). Therefore, the sum of its residues at the punctures is 
equal to the sum of other residues with negative sign. There are poles of two types. First 
of all, Q has poles at the poles 7 S of C. Note that 5C has pole of the second order at 7 S . 
Taking into account that C* has zero at 7 S we obtain 

res 7s fi =< C*5LC > A5z(-f s ) = 5k(^ s ) A 5z(j s ). (3.5) 

The last equality follows from the standard formula for a variation of the eigenvalue of an 
operator. The second term in ( ^.1| ) has the same residue at 7 S . 

The second set of poles of Q is a set of the branching points qi of the cover. The pole of 
C* at qi cancels with the zero of the differential dz, dz(qi) = 0, considered as differential on 
T. The function C is holomorphic at qi. If we take an expansion of C in the local coordinate 
[z — z(qi)) 1 / 2 (in general position when the branching point is simple) and consider its 
variation we get that 

dC 

SC=-^5z( qi )+0(l). (3.6) 
dz 

Therefore, SC has simple pole at qi. In the similar way we obtain 

dk 

8k = -—Sz(q t ). (3.7) 
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Equalities ( pT6| ) and ( |3.7| ) imply that 



res q U = res qi 



< C*5LdC > A 



5kdz 
dk 



(3.8) 



Due to skew-symmetry of the wedge product we we may replace 5L in ( |3.8| ) by (5L — 5k). 
Then using identities C*(5L — 5k) = 5C*(k — L) and (k — L)dC = (dL — dk)C we obtain 



res gi fi 



-res 9i < 5C*C > A5kdz = res 9i < C*5C > A5kdz. 



(3.9) 



(Note, that the term with dL does not contributes to the residue, because dL{qi) = 0.) The 
right hand side of (|3.9|) cancels with a residue of the second term in the sum (|3.1|) . The 
theorem is proved. 

Now let us express to in terms of the coordinates ( |2.39|) on M.. Using the gauge trans- 
formation 

L = GLG-\ C = GC (3.10) 



where G is given by ( |2.18|) , we obtain 



res 2=0 Tr [5L A 5h + C' 1 (5L A 5C + [6h, L] A 5C - 5h C A 5k 



(3.11) 



where 5h = 5GG 1 , 5h = diag (8xiC(z)). From the relation LC = Ck it follows that: 
Tr (C- I [5h, 1} A 5C) = Tr (p^Sh A {L5C - 5Ck)) = -Tr (p^Sh A {5LC - C5k) 



Therefore, 

u = res 2=0 Tr (25L A 5h + C~ X 5L A 5C) 
The first term equals J2i A 5pi. The last term equals 



(3.12) 



~ Tr Co" 1 ^^) A 5C = ^ Tr (c^iSBjA,) A <5C + 5C, 1 A {B+6Aj)C ) , (3.13) 



where Co is the matrix cf of leading coefficients of the expansions Q2.45| ). From (2.46) it 
follows that 

N N 

£ (A ia 5c| + <L4 J>C *) =0, J (c* ki 6B° + 5c* kl Bf) = , 
i=i j'=i 



where c* fc are matrix elements of C . Substitution of the last formulae into ( |3.13j ) completes 
the proof of the following theorem. 



Theorem 3.2 The symplectic form uj given by 1\) equals 

u = ~ E ( S Pi ^tei + H SBf A 5A ia ) . 

i=l V a=l ) 



(3.14) 
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When the symplectic structure u is identified with a standard, it can be directly checked 
that equations ( 2.28| ), ( 2.35 ) are Hamiltonian with respect to uj and with the Hamiltonian 



1 N 1 

H = ^Y,P 2 i ~ 2 J2( b t a j)( b t a i)p( x i - x i)- ( 3 - 15 ) 

i=l i^j 



Nevertheless, we would like to show that the existence of a Hamiltonian for system ( 2.28|) 



( p.35| ) can be proved in the framework of our approach without use of the explicit formula 
for the symplectic structure. 

By definition a vector field dt on a symplectic manifold is Hamiltonian, if the contraction 
iQ t uj{X) = u(X,dt) of the symplectic form is an exact one-form dH(X). The function H is 
the Hamiltonian corresponding to the vector field d t . The equations d t L = [M,L], d t k = 0, 
and the equation 

d t C{t, P) = MC{t, P) + fjL(P, t)C(t, P), (3.16) 
where /i(-P, t) is a scalar function, imply 



N 



i 9tUJ = ^res Pi (< C*(5L + 6k(z))(M + y)C{z,k) > - < C*(z,k)[M, L]5C >)dz, (3.17) 

i=l 

Note, that if a matrix A (z) is holomorphic outside Pj, then the differential < C*AC > dz 
is holomorphic outside Pj, as well. Therefore, the sum of its residues at Pj is equal to zero. 
Using that and the relations 

< C*[M, L]5C > = < C*M{L - k)5C > = < C*{z, k)M(Sk - 5L)C >, 

we get 

N N 

i Qt uj = ^res Pi (< C*(5L + 5k)C > n(P,t))dz = 2^res Pi 5fc n(P,t)dz. (3.18) 

i=l i=l 

The singular part of the function fi(P,t) is equal to the singular part of the eigenvalues of 
the second Lax operator 

H S {z) = -z~ 2 + z~ 1 2(k j (z)) + 0(1), (3.19) 
where kj(z) is the expansion of k(z) at Pj (see (4.8) in ||15|| ). Hence, 

i dt uo = 2 res z=0 Tr [z~ 2 5k + z^Sk 2 ) dz = 2 <5Tr k 2 = 25 Tr L 2 = 5H. (3.20) 

Now, let us define a symplectic structure for spin generalization of the elliptic RS system 
by the formula 

N 

u = J>es p (< C*(z, k){8L{z)L~ l + 5 In k) A 8C(z) >) dz, (3.21) 

i=l 



where L and C are the Lax operator ([2.56) and its eigenvector. 
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Theorem 3.3 The two-form (\3. 2\ ) equals 



uj = 2 J2 5z( Cis) A 51nfc(7 s ). 

s=l 



(3.22) 



Equation ( \2.5(\ )- fl2.52\ ) are Hamiltonian with respect to this symplectic structure with the 
Hamiltonian 

N 

H = Y / (bfa l ). (3.23) 

8=1 



Note, that ( |3.22|) implies that uj is closed and does define a symplectic structure on M.. The 
proof of Theorem 3.3 goes almost identically to the previous case. 

At this moment, we do not know for / > 1 an explicit expression for uj in the original 
coordinates on j\A. Formula ( [3.21|) contains the inverse matrix L^ 1 . For I = 1 it can be 
written explicitly. Namely, let L be the matrix 



a(z + Xj- Xj) 
1 a(xi - xj - rj) 



which is (up to gauge transformation (|2.18|) and a scalar factor) equal to (|2.56|) for I — 1. 
Then the entries of the inverse matrix equal: 



fr 



_ x a(z + xj — x m — (N — 2)77) Y\ k -t m (y(xj -x k + i]) Y\ k a(xj -x k -rj) 



a(z + r])a(z - (N - 1)77) Eh^ <r(,Xj - x k ) Uk^m &(xm - x h ) 
This formula has been used for the proof of the following theorem. 



Then the form uj defined by ( \3. 21\ ) equals 

uj = 5 In fi A 5xi + V(xi — Xj)5xi A 5x 



Theorem 3.4 Let L be the matrix Lij = /,,$( ), where $ is given by ( \2.55\ ). 

(3.24) 
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'■j 



where V(x) is defined in ( $.4^ )- 



Our approach is evenly applicable to the nested Bethe ansatz equations. It gives Hamiltonian 
version of the proof that discrete evolution a;™ is a canonical transform of the RS symplectic 



structure |3.24j) . This result was obtained for the first time in [|12|] with the help of Lagrangian 
interpretation of (|1.1|). 
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